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We extend the description of gravitational waves emitted by binary black holes during the final 
stages of inspiral and merger by introducing in the third post-Newtonian (3PN) effective-one-body 
(EOB) templates seven new “flexibility” parameters that affect the two-body dynamics and gravita¬ 
tional radiation emission. The plausible ranges of these flexibility parameters, notably the parameter 
characterising the fourth post-Newtonian effects in the dynamics, are estimated. Using these esti¬ 
mates, we show that the currently available standard 3PN bank of EOB templates does “span” the 
space of signals opened up by all the flexibility parameters, in that their maximized mutual overlaps 
are larger than 96.5%. This confirms the effectualness of 3PN EOB templates for the detection 
of binary black holes in gravitational-wave data from interferometric detectors. The possibility 
to drastically reduce the number of EOB templates using a few “universal” phasing functions is 
suggested. 

PACS numbers: 04.3.ODb, 04.25.Nx, 04.80.Nn, 95.55.Ym 


I. INTRODUCTION 

Current theoretical understanding, stemming both from general relativity and astrophysics, places black hole bi¬ 
naries at the top of the list of candidate sources for the interferometric gravitational-wave detectors that are nearing 
the completion of their construction phase. On the one hand, black hole binaries are by far sources whose dynamics 
(early inspiral and late time quasi-normal mode ringing to a large extent, late inspiral, plunge and merger to a lesser 
extent) is better understood than other sources, such as supernovae or relativistic instabilities in neutron stars, so 
that it is possible to construct reasonably good template waveforms to extract signals out of noise. On the other 
hand, astrophysical rate estimates of black hole binary coalescences, though not known accurately, have a range whose 
upper limit is large enough to expect a few mergers per year within a distance of 150 Mpc 0]. 

The merger phase of binaries consisting of two 15M 0 black holes takes place right in the heart of the LIGO- 
VIRGO-GEO sensitivity band giving us the best possible picture of this highly non-linear evolution. Eventually, 
when detectors reach good sensitivity levels, one hopes to learn experimentally about this strong gravity regime which 
has been a subject of intense analytical and numerical studies for more than a decade. In the meantime, what is 
needed is a set of model waveforms or templates that describe the dynamics close to the merger phase accurately 
enough so that only a small fraction (< 10%) of all events will go undetected. 

Two general viewpoints are possible to reach this goal. A maximalist one or a minimalist one. The maximalist 
viewpoint consists in enlarging as much as is conceivable (using in a “democratic” way all available methods in the 
literature for treating binary coalescence) the bank of filters, with the hope that even the methods which appear a 
priori less reliable than others might, by accident, happen to describe a good approximation to the “real” signal. This 
is the viewpoint taken by Buonanno, Chen and Vallisneri in their careful, and detailed analysis in Ref. |2(, on the 
basis of which they advocate expanding the net by using a multiparameter template family able to approximate most 
of the results of the conceivable analytical methods. However, this “democratic” attitude comes at a cost that calls 
for an alternate strategy that we explore here. The problem with this method is that it leads to a dramatic increase in 
the total number of templates from ~ 25 templates to ~ 10 4 templates; which has the bad consequence that it leads 
a larger false alarm rate. (These estimates of the number of templates are those obtained, as in !(2], by “dividing” the 
parameter space by the local span of the template at the minimal match. This is an underestimate because it neglects 
boundary effects. For instance, a more realistic estimate of the number of the third post-Newtonian effective-one-body 
templates would be ~ 150 rather than 25.) 
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By contrast, we advocate here a minimalist viewpoint consisting in (i) focussing exclusively on the best available an¬ 
alytical description, and (ii) generalizing this description by adding several parameters that describe “new directions” 
corresponding to physical effects not perfectly modelled by this description. The most important of these directions 
are the effects due to higher post-Newtonian (PN) effects, not yet calculated, but known to exist. A thorough study 
of the robustness of our preferred description against the inclusion of currently unknown effects should allow an in¬ 
formed and judicious covering of the parameter space of interest without overtly expanding the size of the total bank 
of templates. 

The best available analytical description at present is, in our opinion, the “effective-one-body” (EOB) approach 
proposed and constructed for non-spinning bodies at second post-Newtonian (2PN) order by Buonanno and Damour 
% |I], extended to third post-Newtonian (3PN) order by Damour, Jaranowski and Schafer [|] and generalized to 
spinning bodies by Damour §J]. On the one hand, the EOB maps, using the Hamilton-Jacobi formalism, the real 
“conservative” (in the absence of radiation reaction) dynamics of two bodies with masses m\ and m 2 into an EOB 
problem of a test particle of mass /i = Mr] (where M = m± + m 2 and r] = m\m 2 /M 2 ), moving (essentially) in 
an effective background metric g®® which is a deformation of the Schwarzschild metric with deformation parameter 
T], Further, by supplementing the above dynamics by an additional radiation-reaction force obtained from a Pade 
resummation of the gravitational-wave flux, it allows for the first time the possibility to go beyond the adiabatic 
approximation and to analytically discuss the transition from inspiral to plunge and the subsequent match to merger 
and ringing. The implications of the EOB templates for data analysis of binary black holes were explored in [?] where 
it was shown that the signal-to-noise ratio (SNR) is significantly enhanced relative to the usual PN templates due to 
inclusion of the plunge signal. The EOB formalism does also provide initial dynamical data (position and momenta) 
for two black holes at the beginning of the plunge to be used in numerical relativity to construct gravitational data 
like metric and its time derivative and evolve Einstein’s full equations through the merger phase. 

The analytical prediction of the EOB method (including spin) for invariant functions were compared to numerical 
results based on the helical Killing vector approach Jj for circular orbits of corotating black holes by Damour, 
Gourgouhlon and Grandclement |s] and shown to agree remarkably well. The agreement was robust against choices 
of resummation of the EOB potential and improved with the PN order. Recently, Buonanno, Chen and Vallisneri [§] 
made a detailed and exhaustive comparison of all currently available waveforms for non-spinning binary black holes 
resulting from different approximations. This study showed (among other results) that EOB models are more reliable 
and robust than other non-adiabatic models. 

Recently Blanchet [lcj has made a comparison of the straightforward PN predictions with || and shown that at 
3PN order they are as close to the numerical results as the resummed approaches. While it is indeed interesting 
to note this closeness of the results derived from one particular non-resummed 3PN function (the “energy function” 
E(u>) for circular orbits) to the numerical results (and to the EOB ones), we still do not see any way yet by which, 
as the Hamiltonian H( r, p) does in the EOB approach, the bare PN results for the “ circular energy function” E{u>) 
can be used to define templates beyond the adiabatic approximation. As Ref. has shown the importance of going 
beyond the adiabatic approximation in describing the smooth transition between the inspiral and the plunge, and as 
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has shown how significant was the contribution of the plunge to the SNR, we consider that the EOB waveforms 
are the best, currently available, analytical templates for binary black hole coalescences. 

In this paper, we wish to strain to extremes the flexibility of the EOB formalism by tugging it in directions where 
it can be theoretically pulled and by locating directions it is most likely to yield under deformations by unmodelled 
effects (including higher order PN effects). We shall introduce seven “flexibility” directions i.e. seven new flexibility 
parameters. Then we will investigate the “span” of the original (3PN) EOB templates in the space of waveforms 
opened up by our extension of the EOB waveforms into seven new “flexibility” directions. By span of a given bank of 
templates, we mean the region of signal space which is well modelled by some template in the bank, i.e. the set of signals 
S such that the maximized overlap of S with some template T is larger than 0.965 . Among our seven flexibility 
parameters, some like 65 represent higher order (fourth post-Newtonian, 4PN) corrections in the dynamics, some like 0, 
the arbitrariness in the best available 3PN gravitational-wave flux calculation due to incompleteness of the Hadamard 
partie finie regularization, and others like cp denote a parameter used to factor the gravitational-wave flux in m and 
accelerate the convergence of the Pade approximants to the numerical flux in the test-mass case. Further, the current 
development of EOB has made, at several stages, specific choices of representation of various physical effects and, as 
in any analytical construction, the choices were the simplest that one could apparently make. We then explored the 


1 When we don’t have a perfect template to capture a signal then the effective distance up to which a detector could have ideally seen 
goes down. Suppose with the use of the correct template an antenna could detect sources at a distance D. If the best overlap we can 
achieve with the true signal is / then that distance drops to fD and the new rate of events would be proportional to (/D) 3 . In other 
words, the fractional decrease in the number of events is (1 — / 3 ). By demanding that (1 — / 3 ) < 0.1 (i.e. a loss of no more than 10% 
of all potential events) we get / ~ 0.965. 
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effects induced by a modification of these simple choices, i.e. the consideration of new versions of EOB, characterised 
by different parameter values reflecting other allowed more complex choices. These comprise the remaining four 
parameters and include z-i (£ 2 ) a parameter appearing in the effective Hamiltonian H e s of the EOB, /NonAdiab, a 
parameter to modify the simplest treatment of non-adiabatic effects in the current version of EOB templates; /NonCirc, 
a parameter to modify the simplest treatment of non-circular effects in the current version of EOB templates and finally 
/transition to allow the possibility that the transition between plunge and ringing may occur at frequencies different 
from that assumed in the simplest EOB model. The seven parameters (b 5 , 0, c Pl C 2 , /NonAdiab, /NonCirc, /transition), 
are referred to as the flexibility parameters. Varying them and testing the change of the physical predictions under 
reasonable variation of these parameters is a way of probing the overall robustness of the EOB framework. What 
we specifically investigate is whether a bank of standard 3PN EOB templates is sufficient to represent all plausibly 
relevant extended family of waveforms generated by these new flexibility parameters. 

In the context of this investigation our waveforms will be parameterized by two sets of parameters: (a) The first set 
consists of the usual intrinsic and extrinsic parameters entering the construction of standard waveforms, such as the 
masses of the component objects and their spins, some reference phase, etc., denoted collectively bvjok, k = 1,... ,K. 
In this work we shall only deal with non-spinning point-particles in the restricted PN approximation |l2j which requires 
K = 4 with pi 2 = 7711,2 denoting the masses of the two bodies, P 3 = t re f a reference time (related to the instant of 
coalescence), and P 4 = 4> 0 the phase of the wave at the reference time, (b) The second set consists of the flexibility 
parameters introduced above. We shall denote these flexibility parameters as 7r a , a = 1,..., A, with {A = 7) 

TTl = 65 , 71" 2=0, 7T 3 = Cp , 7T 4 = ,Z 2 (C 2 ), 

775 = /NonAdiab, TTfj = /NonCirc, 7T 7 — /transition- (1 ■ 1) 

For studies of the span of a bank of templates of the kind we propose to do in this paper, it is helpful to introduce 
the notions of a standard or fiducial template and its associated variant or flexed signal constructed by turning on 
one of the flexibility parameters. The term fiducial template is used to represent a waveform constructed in a certain 
approximation and at a given PN order with a fixed set of values of the unknown parameters introduced above. In this 
paper, our fiducial template will be the standard EOB waveform (see Sec. IIIII) at the 3PN order with the flexibility 
parameters all set to zero: 


Fiducial Template = h(t;pk,^b = 0,V6) = T(t;p k ). (1.2) 

In contrast, the associated flexed signal will again be the EOB waveform at 3PN order with all but one of the flexibility 
parameters 7 r 0 set to zero: 


7 r a -Flexed Signal = h(t;p k ,n a £ 0 , 7 r b = 0,V b ^ a) = S(t;p k , n a ). (1.3) 

In other words, in our test of robustness we do not allow all the seven parameters to vary simultaneously. Such 
a variation would lead to a formidably high dimensional parameter space which is computationally impossible to 
investigate at the moment. Rather, our aim is to study the effect of each flexibility parameter independently and 
to gauge the extent to which our standard fiducial template waveform can mimic the changes brought about by 
the flexibility parameters n a by a mere variation of the intrinsic parameters p k ■ Such a systematic study allows us 
to isolate and identify the most important unknown physical effects, and decide if it is necessary to introduce the 
corresponding flex parameter as an additional parameter in search templates used in the detection of gravitational 
waves from black-hole binaries. 

As in earlier work, our main tool for measuring the span of a bank of templates is the overlap of a standard fiducial 
template with a given flexed signal. (See Sec IVTI for the definition of the overlap.) In this study, we use two measures 
of “good overlaps”: faithfulness and effectualness 0 - A template is said to be faithful if its overlap with a flexed signal 
waveform of exactly the same intrinsic parameter values is larger than 0.965 (after maximization over the extrinsic 
parameters). It is expected that faithful templates are also good at estimating source parameters although this is not 
guaranteed to be the case. A template is said to be effectual if its overlap with a flexed signal waveform, maximized 
over all the intrinsic and extrinsic parameters , is larger than 0.965. Obviously, every faithful template is necessarily 
effectual but not all effectual templates are faithful. Note that the notions of faithfulness and effectualness might 
depend on the particular flexibility direction which is explored: While a template waveform could be faithful with 
respect to the 7Ti-flexed signal, it might only be effectual with respect to the 7r 2 -flexed one and neither with respect 
to the 7 T 3 -flexed one. 

Implementing the above analysis we conclude that the standard 3PN EOB templates are effectual with respect to all 
flexibility parameters introduced in this study including the parameter 65 characterising the 4PN dynamical effects. 
In other words, the span of the bank of 3PN EOB templates is large enough to cover the space of signals described 
by the physically plausible ranges of the seven flexibility parameters considered here. No additional extra parameters 
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are required in the detection templates to model the more complex choices possible in the EOB approach at the 3PN 
level or the dominant dynamical effects at the 4PN order. In particular, there is no need to increase the total number 
of templates beyond the level required for the standard 3PN EOB templates. 


II. THIRD POST-NEWTONIAN DYNAMICS AND ENERGY FLUX 


The conservative dynamics of binary systems in the PN approach has now been determined to 3PN accuracy. Two 
independent calculations, one based on the canonical ADM approach together with the standard Hadamard partie 
finie regularization for the self-field effects mill and the second, a direct 3PN iteration of the equations of motion 
in harmonic coordinates supplemented by an extended Hadamard partie finie regularization jlj, I16l U7I agree that 
the 3PN dynamics and consequently conserved quantities like energy are fully determined except for one arbitrary 
parameter called u> s in the ADM approach and A in the harmonic coordinates related by, 


3 1987 

_ TT Ws ~~ 3080' 


( 2 . 1 ) 


The Hadamard regularization of the self-field of point particles used in 0 0 E13 E11 has the serious 
drawback of violating the gauge symmetry of perturbative general relativity (diffeomorphism invariance), and thereby 
of breaking the crucial link between Bianchi identities and equations of motion. This explains why the Hadamard- 
based works in iii m ee in m were unable to fix the parameter u> s . Recently, Damour, Jaranowski and Schafer 
ESI have proposed to use a better regularization scheme, one which respects the gauge symmetry of perturbative 
general relativity: dimensional regularization. They have implemented this improved regularization scheme, which 
led them to a unique determination EE] of the parameter u> 8 , namely uj s = 0, corresponding to A = —1987/3080. 
Thus the conservative dynamics is completely determined to 3PN order within the ADM approach using dimensional 
regularization. Though it will be interesting to reconfirm the value of A by other treatments, we believe that the 
result of Hi is trustable especially in view of the obtention there, by the same regularization method, of the unique 
Poincare-invariant momentum-dependent part of the Hamiltonian. Thus, for all applications including data analysis, 
there is no arbitrariness in 3PN dynamics, and consistent with this, in this paper we set io s = 0 or equivalently 
A = -1987/3080. 

On the other hand, the gravitational-wave energy flux from binary systems has been computed using the multipolar 
post-Minkowskian approach .2Qj in harmonic coordinates and Hadamard regularization to 3.5PN accuracy. Unlike 
at earlier orders l2ll| the instantaneous and hereditary contributions do not remain isolated. At 3PN order, in 
addition to the instantaneous terms, the tails-of-tails and tail-squared terms also contribute. Fortunately, they have 
been computed by Blanchet EE] who has also computed the tail contribution at 3.5PN order. The gravitational 
wave energy flux contains the 3PN-accurate time derivative of the mass quadrupole moment leading to a specific 
A dependence which as explained earlier is now known since cu s is computed. However, the incompleteness of the 
Hadamard regularization introduces additional arbitrary parameters in the mass quadrupole moment leading to three 
new undetermined parameters that combined into the unique quantity 9 in the circular energy flux. Unfortunately, up 
to now no alternate regularization or calculations without regularizations exist that provide the value of 9. Thus, one 
has to reckon with this arbitrary parameter in the templates that one constructs and the best one can do is estimate 
its implications for data analysis of inspiraling compact binaries as in [3j or in the present work. 

The expression for the 3PN energy function El El El and 3.5PN flux function EE |24j and the resulting 3PN 
and 3.5PN coefficients in various phasing formulas discussed in Jjj are summarised in 1 2 "I . Though they are the basis 
of the present analysis, they are not reproduced here for reasons of brevity and we refer the reader to l2f| for those 
expressions. 

In this work, we have used the 3PN-accurate flux function because the standard near-diagonal Pade of the 3.5PN 
flux function involves a spurious pole in the physically relevant range of variation of v. We leave to future work an 
investigation of alternative Pade’s of the 3.5PN flux free of such spurious poles. In view of the numerical smallness of 
the 3.5PN contribution to the flux, we expect no significant change in our physical conclusions. 


III. TRANSITION FROM INSPIRAL TO MERGER—THE 3PN EFFECTIVE-ONE-BODY MODEL 

The starting motivation of the EOB approach is to try to capture in a small number of numerical coefficients the 
essential invariant PN contributions from among the plethora of terms that exist in the complete PN expansion of 
the binary’s equations of motion, in the belief that many of these terms are gauge artefacts and hence irrelevant. It 
is also strongly motivated by the need to look for an analytic route to go beyond the adiabatic approximation which 
breaks down before the last stable orbit. We recall that the standard PN treatments based on invariant functions 
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(. E(w ), are limited by the adiabatic approximation (and cannot describe the transition to plunge), while the 

treatments based on the direct use of (non-resummed) PN-expanded equations of motion are unreliable [3 because 
of poor convergence of the straightforward PN-expanded equations of motion. 

As shown in |^, at 2PN order the mapping to EOB is eventually unique (when imposing some general requirement). 
The waveform, the equations governing the evolution of the orbital phase and the initial conditions to integrate them 
through the plunge are discussed in H) and were used in 0 to construct 2PN EOB templates and investigate their 
performance. At 3PN order, on the other hand, the situation is more involved. When requiring that the relative 
motion be equivalent to geodesic motion in some effective metric, there are more constraints than free parameters in 
the energy map and effective metric. This led la to an extension of the 2PN EOB construction (non-geodesic motion) 
involving a larger variety of choices. In Ref. [5] the following generalized 3PN EOB Hamiltonian was introduced: 


Heff(r,p) = \ A(r) 


A i r ) 

Dir) 


- 1 ) (n ■ p ) 2 + —( 2-1 (p 2 ) 2 + z-i p 2 (n ■ p ) 2 + 23 (n- p) 4 ) 


(3.1) 


where the functions A{r) and D{r) are given by the components of the effective spherically symmetric metric <?®®: 
A{r) = — ffoo( r ) an d D{r)/A{r) = ff®®(r) (they also depend on the parameters z\ and 22 ; see below). Here r and p 
denote the (scaled) canonical coordinates of the effective dynamics, r = |r|, 11 = r/r; r is dimensionless, being scaled 
by GM. The effective position vector r is linked id to the relative position vector xi — X 2 of the two holes in ADM 
coordinates by a post-Newtonian expansion which starts as: r = (xi — X 2 )/{GM) + (P(c _2 ). 

The parameters 21 , 22 , and 23 are arbitrary, but subject to the constraint 


821 + 422 + 323 = 6(4 - 377 ) 77 , 


(3.2) 


which forbids the geodesic choice 21 = 22 = 23 = 0 of the 2PN EOB, but allows the minimally non-geodesic choice, 
21 = 22 = 0 and 23 = 2(4 — 377 ) 77 . [See, however, below our discussion of non-minimal choices as flexibility directions]. 

In spherical coordinates (r, 9, </>), restricting the motion to the equatorial plane 9 = 7 t/2 , the Hamiltonian Eq. 113.11) 
can be written as 


H eS (r,p r ,p ! f,;zi,Z2,z 3 ) = 


\ 


4 


A(r; 21 ) < 1 + A(r; 2 i)Z)(r; 21 , 2 2 ) 1 P 2 + + z (r,Pr,P<f>', Zi, 22 , 23 ) > , (3.3) 


where 


Z{r,Pr,p ( j,]Zl,Z2,Z 3 ) = 


Zl [Pr + ^JJ + z 2 Pr+ Z 3Pt 


The functions A(r\z\) and D(r\ Z 2 ) 1 depend on z\ and 


A(r;zi) = 1 
D(r; 2 i, 2 2 ) _1 = 1 


2 277 04(21) 


677 2(26 — 377)77 — 72i ^ 22 

r 2 ^ r 3 


where 


. , .'94 41 2 

04 ( 21 ) = ( y - ) V - 21 . 


(3.4) 


(3.5a) 

(3.5b) 


(3.6) 


The (scaled) 3PN EOB-improved real Hamiltonian is the following function of the EOB Hamiltonian Eq. (13.31) : 


-ffreal = “ \] 1 + 2?7 {H e ff — 1) . 
77 


(3.7) 
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The equations of motion have the form of the usual Hamilton equations: 


dr 

dt 

dcj) 

dt 

dp r 

dt 

dptf, 

dt 

where is the </> component of the damping force. 



(3.8a) 

dp r 

^^real 

dp$ 

(3.8b) 


(3.8c) 

dr ’ 

F<t> > 

(3.8d) 


A. Pade approximants of A 


The straightforward PN expansion of the function A, in terms of the variable u = 1/r, reads: 

A(u) = 1 — 2u + 2rju 3 + ai(if)u A + 05 {rj)u 5 + 0(u 6 ). (3.9) 

To improve the convergence of the PN expansion ESJ we introduce the following sequence of Pade approximants of 

A0: 


Aipn(m) = 

1 — 2 u , 

(3.10a) 

A2Pn(«) = 

1 - ( 2 - h ) u 

(3.10b) 

1 + ^7777 + 77U 2 ’ 

A3PN(0) = 

2(4 - 77) + (04(77) - 16 + 877)77 

(3.10c) 

2(4 - 77) + (04(77) + 477)77 + 2(04(77) + 477)zt 2 + 4(04(77) + 77 2 )t7 3 ’ 

A4PN (0; 05) = 

(16 — 877 — 04(77)) — (32 — 2477 — 404(77) — 05(t7))t7 

(3.10d) 

d4PN(o; 05) 

where 



d 4PN (0; o 5 ) = (I6-877- 

04 (77)) + (877 + 204(77) + 05(77)) (77 + 2t7 2 ) + 2(8?7 2 + (4 + 77)04(77) + 205(77)) 

77 3 

+ (16t7 2 + ai ( rj ) 2 + 87704(77) + (8 — 277)05(77)) u 4 . 

(3.11) 


Reference 0 has studied some variants of the specific Pade choices made in Eqs. emu and found that they had very 
little effect on physical quantities down to the last stable orbit. Therefore, we shall not include below, among our 
flexibility directions, the ones corresponding to different choices of definition of A function. We shall always define it 
using the P „-type Pade used above. While the 3PN-level coefficient 04 ( 77 ) = 6477 is known when using w s = 0 0 [see 
Eq. iTHbll below], the 4PN-level coefficient 05 ( 77 ) introduced here is unknown. Its possible values will be discussed 
below. 


B. 3PN EOB adiabatic initial data 

The initial dimensionless frequency u)q depends on the initial frequency /g ff of the gravitational wave and the total 
mass M of the binary system: 


~ _ GM7T/g w 
W °- -?- 


(3.12) 
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In the following equations A and D are treated as functions of it, A' = dA/du. The initial value of ro = 1/ito one 
obtains solving numerically equation: 


So = u 3 ' 2 


-U' 


l + 2r7 


A 


\ \ \J A + \uA' 


- 1 


The initial momenta are then obtained from equations: 


P° = 


Pr = 


A! 


y it (2 A + uA') 

u (2 A + uA') A'D 
A [2u(A') 2 + AA! - uAA"} 


o) 


cu 0 


(3.13) 


(3.14a) 

(3.14b) 


C. 3PN EOB light ring 

The light ring coordinate iiiight ring = 1/flight ring is the solution of equation 


A(u) + -uA\u) = 0 . 


(3.15) 


See Ref. 0] for a discussion of the physics associated with the light ring within the EOB framework. [See also [23 for 
alternative views based on the non-resummed 3PN energy function E(ui).] 


IV. ESTIMATING THE EFFECT OF UNKNOWN PHYSICS 

Our main goal is to vary the flexibility parameters within a certain range motivated by physical arguments to be 
discussed in Secs. ITvTniTvTn and to determine the degradation caused by such a variation on the overlap of flexed 
waveforms with fiducial waveforms. If the degradation is small then we further explore the maximum extent to which 
the parameters can be meaningfully varied so that the effectualness (see Sec. I VI Al for a definition) still remains more 
than 96.5%. We limit the range of variation of each parameter so that the effective potential, binding energy and 
energy flux remain regular and meaningful for values of the parameter in that range. The natural range over which 
the parameters are expected to vary is summarized in the first row of Table rm 


A. Higher order PN dynamics parameter b$ 


In terms of the inverse u = 1/r of the radial effective coordinate r the PN expansion 13.911 of the function A(u) = 
~ffoo( u ) can be written as: 


A(u) = 1 — 2 u + 63 77 u 3 + 64 77 u A + 6577 [1 + 0(r])\ u 5 + 0(u 6 ), 


(4.1) 


where 


— 2, 


(4.2a) 


&4 


94 41 , 

-7T 

3 32 


(4.2b) 


We have included in Eq. JUJ the information that for the 2PN (oc u 3 ) and 3PN (oc u 4 ) levels there has been rather 
miraculous cancellations to leave only terms linear in r). Moreover, we also know that to all orders (starting from 
2PN) the terms oc rf vanish. We expect that the terms linear in 77 in the higher PN coefficients dominate over the 
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nonlinear ones. In particular, we expect that in the 4PN coefficient 05 ( 77 ) = 6577 + 0(rj 2 ), the term 0(rj 2 ) can be 
neglected. Finally, we shall work under the simple assumption < 25 ( 77 ) = 6577 , i.e. with 

A(u) = 1 — 2 u + 63 77 u 3 + 64 77 u 4 + 65 rju 5 + 0(u 6 ). (4.3) 

The main aim of this subsection will be to estimate the plausible order of magnitude of the 4PN coefficient 65 . 

We can start to guess a plausible range of values of 65 by the following reasoning. It is plausible to expect that the 
function A(u) be a meromorphic function of u (or, at least, be close to a meromorphic function). The growth with n 
of the Taylor coefficients b n of a meromorphic function is determined by the location of the nearest singularity in the 
complex plane of the function A(u). If the nearest singularity is located at u = 1/6, the Taylor coefficients b n of A(u) 
behave, when n increases, roughly proportionally to b n . We can always parametrise this behaviour (without loss of 
generality) as 


b n ~ k b n ~ 3 . 

Using Eq. (TT~fl we can deduce b and k from £>3 and 64 : 

k ~ b 3 = 2 , 
b ~ 54/63 — 9.3. 


(4.4) 


(4.5a) 

(4.5b) 


This yields the guess 


65 ^ 170. (4.6) 

Note that the values of b and k from Eqs. 14.511 gives also a “prediction” for 62 which is 

b 2 ~ k/b ~ 0.2. (4.7) 

This type of value is small enough not to make a physical difference from the exact value b 2 = 0, and moreover it is 
compatible with the fact that the detailed calculation of b 2 gives a cancellation of the type b 2 = 0.25 — 0.25 (see text 
below and Tabled, each term being indeed of order 0.2. This is consistent with a power-law growth of the typical 
contributions entering b n . 

A first guess is therefore that 65 is positive (because this is true for 63 and 64 ) and smaller than 200 (in round 
numbers). To go beyond this guess we studied in detail the various PN contributions to the successive b n , with the 
aim of detecting a pattern. We explain in detail our study in the remainder of the subsection. 


1. E = E(j) for circular orbits 


We work here with the Hamiltonian describing the real relative motion of the two bodies in their center-of-mass 
reference frame (the superscript NR denotes a “non-relativistic” Hamiltonian, i.e. the Hamiltonian without the rest- 
mass contribution). It reads 


1 


1 


1 


1 


# (r, p) = Hn (r, p) + -jH 1PN (r, p) + -tH 2PN (r, p) + -^H 3PN (r. p) + H iPN (r, p) 


(4.8) 


where 


xi - x 2 

r =- 

GM 


P 


Pi 


P2 j/NR _ -^ NR 
M _ M 


(4.9) 


Let us note that in this subsection (and only here) r and p denote the canonical coordinates of the real (relative) 
two-body dynamics. 

Circular motion is defined through the condition 


n ■ P = 0, 


(4.10) 
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where n = r/r and r = |rj. Under the condition Eq. iITTol) the Hamiltonians -Bn through i? 4 PN have the structure 


ff _ P 2 1 

Hn - ~2 r ’ 

(4.11a) 

»** - 

(4.11b) 

k—0 


^ 3 (n 2 \ k 

H 2 pn = Y^h 2 r3 _ k , 

(4.11c) 

k—0 


_ 4 (r, 2 \ k 

H 3 PN = J2 h 3 r 4 -k > 

(4.lid) 

k—0 


5 (n 2 \ k 

H *PN = r5 - k ■ 

(4.lie) 

k=0 


The momentum squared p 2 can always be decomposed as 


j 2 

P = (n p) + (n x p) 2 = p 2 r + 

(4.12) 

where p r = n ■ p and 



(4.13) 


Here J is the conserved total angular momentum of the binary system in the center-of-mass reference frame. 

The Hamiltonian iTTSl . after replacing p 2 by j 2 /r 2 , becomes a function of j and r only. This function has to fulfill, 
by virtue of the equations of motion (for circular motion p r = 0 ), the condition 


dH NR (j,r) 

dr 


= 0 . 


(4.14) 


Equation 114.141 gives the link between r and j along circular orbits. We have iteratively solved Eq. 631 for r as a 
function of j, and have substituted this into the Hamiltonian 63- We thus have obtained the relation, valid along 
circular orbits, between the center-of-mass energy E = 1J NR (r, p) of the system and the system’s angular momentum 
j. To simplify displaying this relation we show it with the coefficients h\ of the 1PN Hamiltonian Hipn replaced by 
their explicit general relativistic values. Then the formula reads 


E = — 


2f 


X +4 (9 ' 


V)r 2 

3 


' 3 N 

16-2 £>2 
v k=0 j 


152 — 24 77 — 8 j> + 3)&§-2£>* 


k —0 


k =0 


+ 1700 — 584 rj + 36p 2 + 4^^(fc + 3)(3 77 — 27 — 4 fc) h\ + ( ^(fc + 3)/i£ 


k —0 


\k —0 


- 8 ^(fc + 4 )h\ 
k —0 




(4.15) 


2. EOB potential A(u) calculated from E = E(j) 


In the effective-one-body approach the real “non-relativistic” energy E is the following function of the effective- 
one-body radial potential Wj (u): 



(4.16) 
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where 


Wj(u) = A{u) (l+jV) . (4.17) 

The function A(u ) has a perturbative expansion in u: 

A(u) = 1 + di U + 02 u 2 + 03 u 3 + 04 u 4 + 05 u 5 . (4-18) 

Along circular orbits the effective radial potential Wj(u) attains its minimal value, 

^- Wj (u)=0. (4.19) 

We have iteratively solved Eq. 14.1 91) for u as a function of the small parameter l/j 2 and have substituted this relation 
into the right-hand side of formula tTTTCll . Next we have again expanded the right-hand side of (TCTl in l/j 2 . In such 
a way we have obtained the relation E = E{j) predicted by the general EOB function 14.181) in which the coefficients 
at different powers of l/j 2 depend on the numbers a n entering the function A(u). By comparing these coefficients with 
the respective coefficients of the expansion (IT45l> we are able to (iteratively) express a n in terms of the coefficients 

h^ of the Hamiltonian (TOl) . _ 

After this matching between the generic Hamiltonian iTTSli and the guessed EOB expression eng, each of the 
numbers a n can be represented as a sum of terms which depend on the coefficients of the different PN Hamiltonians. 
E.g., a 2 = a 20 + 021 , where 020 depends only on the coefficients of the Newtonian Hamiltonian and 021 depends 
on the coefficients of the Newtonian and the 1PN Hipn Hamiltonians. More generally, we have 

n—1 

a n = ^ ^ C^nk ? (4.20) 

k —0 


where a n q depends only on H n, a n i depends on H n and H ipn, a n 2 depends on ^ipn an d i^ 2 PN? etc. The 
values of the different coefficients a n k are as follows (here again, to simplify formulas, the coefficients h\ of the 1PN 
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Hamiltonian H ipn 
OlO 

020 

021 

030 

031 

032 

040 

041 

042 

043 

050 

«51 

«52 

053 

054 


have been replaced by their general relativistic values) 


- 2 , 

(4.21a) 

\ (9 + r?) , 

(4.21b) 

(9 + 77 ) , 

(4.21c) 

^(-27-12W) , 

(4.21d) 

^ (-67 + 30 77 + r? 2 ) , 

(4.21e) 

2 X> 2 , 

fc— 0 

(4.21f) 

^ (54 + 72 77 — 11 t ? 2 + 2 77 s ) , 

(4.21g) 

Ttt (-1973 + 1301 77 — 55 r 7 2 77 s ) , 

64 

(4.21h) 

3 

5^(3 — 77 + 8 / 3 ) , 

fc =0 

(4.21i) 

2^3, 

(4-21J) 

1 ( 243 1080 77 + 210 ? 7 2 104t? 3 + 21 t? 4 ) , 

(4.21k) 

1 ( 201017+ 250392 77 46630 rj 2 + 536 t ? 3 + 15 ? 7 4 ) , 

(4.211) 

, 3 3 3 

- ( 207 + 384 k + 128 k2 ~ ( 213 + 128 k) 77 + 7 ? 7 2 ) h\ - ^ k ^(6 + k) h\ , 

° k= 0 fc=l fc =0 

(4.21m) 

4 

5^(5 — 77 + 8 / 3 ) 

(4.21n) 

2 1>4- 

(4.21o) 


k -o 


Because the Hamiltonians from Newtonian through 3PN are completely known, the coefficients ai through <24 are 
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TABLE I: Numerical values of the coefficients b n k- 


n 

bn 

bnO 

bnl 

bn2 

bn3 

bn4 

2 

0 

0.25 

-0.25 




3 

2 

-0.75 

1.875 

0.875 



4 

18.6879 

1.125 

20.3281 

-17.125 

14.3598 


5 

? 

-1.05469 

244.523 

-324.82 

194.214 

? 


fully known. They read 


ai 

= -2, 

(4.22a) 

a 2 

= 0 , 

(4.22b) 

03 

= 277 , 

(4.22c) 


/94 41 


04 

= y-srk 

(4.22d) 


Let us note that 03 and <24 are both proportional to 77 . Many remarkable cancellations occurred to cancel the terms 
proportional to rf and 77 . As for the 4PN-level coefficient 05 = its first three partial contributions < 251 , 052 , 

<153 are known, but its last coefficient <254 is unknown. <25 is a polynomial of order at most 4 in 77 with vanishing term 
oc 77 ° (this can be explicitly checked as the two-body Hamiltonian in the test-mass limit 77 = 0 is known up to all 
orders). As we said above, we expect, as it is the case at lower orders, that in 05 the term oc 77 1 will dominate. Let 
us denote 


a nk =KkV + C>( 7 ? 2 ) j 72 > 2 . (4-23) 

The numerical values of the parameters b n k are given in Tabled] After studying the various possible patterns exhibited 
by Table d we decided to focus on the fact that the column 6 nl of Tabled seems to give a good approximation of the 
final, total value of b n . This would suggest 65 ~ 250, as possible value. 

Let us now display a more explicit form (for the parts which are known) of the 4PN coefficient as = X^fc=o a 5 fc- To 
do this we replace in Eqs. 1 IT 2 H g 21 ol) the 2PN and 3PN coefficients and h 3 by their general relativistic values. 
Out of the 4PN coefficients h\ the leading kinetic term h\ it is fully known (as it is given by the expansion of the free 
Hamiltonian ^ f a vpI + tttJ), for the rest of the terms only their parts oc rf are known (they describe the test-mass 
limit of the two-body dynamics). We parametrise our ignorance of the parts 0(rf) by introducing some quantities 
k = 0,..., 4. Thus we can write 


-T 6 +7 7 *°’ 

(4.24a) 

105 , 

^-+77X4, 

(4.24b) 

105 2 

^- + 77X 4 , 

(4.24c) 

13 3 

Y + ^ ^4 . 

(4.24d) 

45 4 

128+7? 

(4.24e) 

-^-(1 - 3 27 )(1 - 677 + 9 t7 2 - 3 77 3 ). 

256 

(4.24f) 
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Collecting all this partial information together one gets 


571 197 2 \ /41 2 1885 

05 = i —--&* hMer 


\ 2 27 o 35 4 v~^ 

W -W + Ta t i + 2? ?E 


96 / 16 ' 64 


Xi 


fe=o 


(4.25a) 


= 112.3701 ?7 — 13.3127 rj 2 — 1.6875 rf ± 0.5469 rj A ± 2rj \4 ■ 


(4.25b) 


k -0 


Note that the expression confirms that the terms oc rj 2 , rf and rj 4 are sub-dominant. 

Based on the above results we guess the range of plausible values for the parameter 6 5 to be [0,250]. However, 
while exploring robustness we would like to vary 65 beyond this reasonable range subject to the condition that the 
potential remains regular. This condition implies that 65 > —50 as smaller values of 65 introduce poles in the Pade 
approximated version of A(u). Note, however, that all the known successive PN approximations suggest that the 6 n ’s 
are all positive so that the consideration of negative values of 65 test robustness against extreme behaviour of the 
potential. 


B. Location of the pole in energy flux cp 


In Ref. 0 it was argued that we should expect a (simple) pole in the flux function as a function of 


v = (Mcu) 1/3 , 


(4.26) 


where w is the orbital frequency, at the location of the light ring. It was further shown that factoring out the pole 
from the post-Newtonian expansion of the flux before constructing its Pade approximation accelerates convergence to 
the fits to the numerical flux. What happens when we “flex” the position of this pole away from its known test-mass 
value, or its conjectured 77 -dependent 2PN location? In the test-mass approximation, that is 77 —s- 0, the location of the 
pole in the flux function is at Wp ole = 1/ V%. When 77 is different from zero Ref. [nj argued that a good approximation 
to the location of the pole is given by: 


v 


v 

pole 


1 

7 ! 


/ i \ 1/ 2 
1 + 0 7 ? | n 

-gH - ^pole [! + 0.16 (477)] . 

V-mV 


(4.27) 


The location of the pole can significantly change the value of the Pade approximant of the flux function in the 
physically relevant region of the variable v (see below). For this reason it is important to move the pole away from 
its predicted value and assess how such a shift would affect the detectability of the signal. In this work we modify 
the location of the pole by introducing the parameter cp: 


cp _ pole 
pole 1 + C P ' 


(4.28) 


Based on the fact that Vp ole differs, when 77 = 1/4, from the test mass value Pp ole by ~ 16%, an a priori plausible range 
of variation of cp is ±0.2. We also explored larger variations of cp, namely in the range [—0.5, ±0.5], which in the 
comparable mass case amounts to varying the original 2PN pole from 0.6907 in the range [0.4605,1.3814]. Actually, 
values of cp smaller than —0.2 seems to have little effect on the overlaps. [Note that when cp tends to —1 this pushes 
the pole to 77 po ie —> ± 00 .] If cp is taken to be greater than about 0.5 then the location of the pole in the flux will be 
at v < v\ so so that we will not be able to compute the phasing of the waves. This is why we restrict the values of cp 
to be smaller than about 0.5. 


C. Unknown third post-Newtonian energy flux 9 

To estimate the possible range for the unknown parameter 9, let us go back to Ref. [24j where this arbitrariness is 
pointed out and discussed. The parameter 9 is the linear combination 9 = £ ± 2k ± ( of the three coefficients £, k 
and associated with three different kinds of terms. From Eq. (10.6) of Ref. l2~j and more explicitly from work in 
progress ^3 (which uses the generalized Hadamard regularization of Ref. 0 ) it follows that the value of 9 contains 
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a log term in addition to a term of approximate value — 1. This motivates us to suggest that a variation range by a 
factor of order 10 (with both positive and negative signs) is a very generous range for 9 , which is expected to be “of 
order unity”. 

How large can the magnitude of 6 be without introducing any spurious poles in the P-approximant of the flux? 
The answer is that the variation of 6 is bounded from below at 9 = —5 because for 9 < —5 there is a spurious pole. 
However, for values 9 > 0, even as large as 10 5 , there seems to be no irregular behaviour of the P-approximant flux. 
Thus, in our test of robustness we take the minimal range of 9 to be [—5,10] and we also explore the values of 9 > 10. 


D. Modification of the two-body Hamiltonian: z 2 or £2 

The 3PN extension of EOB opened the possibility of introducing two free parameters z\ and z 2 . It is clear that 
taking a non-zero value of z\ goes against the spirit of the EOB resummation, because it takes away a part of the 
basic EOB radial potential A(u) to replace it by a modification of the “centrifugal part” of the potential. [See Eqs. 
inoii inniii above.] Distributing the 3PN effects between A(u) and the centrifugal potential is undesirable because it 
goes against “resumming” all effects in one object: namely A(u). Therefore we continue, as in Q, to fix z\ = 0 and 
this choice simplifies the constraint Eq. ( 13 ^ 1 ) to 


4;z 2 + 3^3 = 24 



V- 


(4.29) 


This leaves us with only one 3PN flexibility parameter linked to this possibility and it is convenient to parametrise it 
by introducing a (2 such that, 

22 = - C 2 z J 3 , 

23 = (1 - C2) 4 , 

where 

4 = 2 (4 - 377) ■ (4.31) 

Making use of z\ = 0, and the above parametrisation, the Hamiltonian Eq. iTHl) reads 


(4.30a) 

(4.30b) 


H e s(r,p r ,p ct> -,C2) = 


\ 


A(r) { 1 + A(r)D(r;C 2 )~ 1 P^ + ^ + : 4 


,2 „2 


/ 1 4 ,3 Pp. 

I 4 Pr + 4 r 2 


(4.32) 


where 


2 2i) 84 /94 41 2 , 

■ 1(r) = 1 ~-r + T* + T‘- = (- 3 - 32 ' I". 


D(r,b)-' = 1 + 4 + 


6p t 4(26 377)77 3(4 — 3t7)C 2 

2 + ^3 


(4.33a) 

(4.33b) 


The values of (2 equal to 0 and 1 correspond to the simple choices z 2 = 0 and 2:3 = 0, respectively. Its variation is of 
order unity to cover this interval and we take its natural range to be [— 2 , 2 ]. 


E. Flexibility parameter for non-adiabaticity /NonAdiab 

The current version of EOB templates chooses the simplest treatment of non-adiabatic effects. In the present study, 
we would like to look at a modification of this choice and to this end we introduce the parameter /NonAdiab hr the 
expression for the angular damping force T,p appearing in Eq. Il3.8dl . More precisely, we modify the current “minimal” 
radiation reaction using: 


> Feb 


1 /NonAdiab 



(Au 2 y 

A' 



(4.34) 
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The combination of factors factored by /NonAdiab vanishes in the adiabatic approximation [see Eq. Il3.14al) l. The 
assumption of adiabaticity is valid for most of the inspiral regime and deviates from it only close to the last stable 
orbit. The modification iTHHli is a simple way of parametrising the effect of different choices in the definition of 
for orbital motions which start deviating from adiabaticity. One generally expects /NonAdiab to be a parameter of 
order unity and it suffices a priori to vary it in the range [—1, +1]. While this is the primary goal, we explore a larger 
range of /NonAdiab in our study of robustness. 


F. Non-circular orbits /NonCirc 


The current version of EOB templates also uses a simplest treatment of non-circular effects, which we would like 
to re-examine here. This is accomplished via another flexibility parameter /NonCirc in the angular damping force. We 
modify the force in a manner similar to the previous case; that is we use 


—> T c/> 


1 T /NonCirc 9 

P<f, U 


(4.35) 


As in the previous case, varying /NonCirc in the range [—1,+1] is expected to be a plausible way of mimicking non- 
minimal choices of the definition of for orbital motions which start deviating from circularity. However, we do 
explore a larger range of /NonCirc in our study of robustness. 


G. Transition between plunge and coalescence /transition 

In the EOB approach the equations representing inspiral are continued through the plunge and eventually matched 
to the quasi-normal modes of the final black hole near the light ring. The exact frequency where this happens cannot 
be decided by the formalism and hence one would like to examine the effect on the waveform by changing the frequency 
of transition between plunge and coalescence and subsequent ringing or, equivalently, stopping the plunge at a point 
different from the point where the EOB waveform is naturally terminated. In the EOB approach the waveform is 
naturally terminated when the radial coordinate gets close to the value rn g ht ring given by a solution to Eq. 13.151) . In 
the current paper we alter the radial location of the light ring by introducing the parameter /transition given by: 

flight ring * flight ring (1 T /transition) ■ (4.36) 

Negative values of the parameter /transition are, rather meaningless since the EOB approximation is expected to 
break down for values of the radial coordinate less than r h g ht ring- We therefore allow only positive values and vary 
/transition in the range [0,1]. Note, however, that the variation in this parameter is going to seriously affect those 
systems which merge in a detector’s sensitivity band since a positive value for this parameter means that we will 
in effect be discarding power in the final phase of the signal. Indeed, in this work, we do not match the plunge 
waveform to the quasi-normal mode expected to ensue soon after. This is because our earlier work in Ref. Q has 
shown that these modes do not contribute significantly to the SNR for those systems whose plunge occurs in the 
detector’s sensitivity band. 


V. HOW ROBUST ARE EOB TEMPLATES? VISUAL COMPARISON 

In this section we discuss the robustness of EOB waveforms by comparing the standard fiducial 3PN EOB templates 
with flexed waveforms constructed by turning on the flexibility parameters discussed in the previous sections. We 
make two different types of comparisons to gauge the extent to which various unknown flexibility parameters at third 
and fourth post-Newtonian orders might affect the dynamics of the two bodies and the radiation they emit. Our 
first comparison consists of a visual inspection of the behaviour of the relevant physical quantity when a particular 
flexibility parameter is varied. This gives us an idea of the nature and the extent of the variation involved while 
testing robustness. Our second comparison goes beyond qualitative tests of robustness by quantitatively measuring 
the span of EOB templates. More precisely, it consists of the computation of the faithfulness and effectualness of the 
fiducial EOB template with the flexed waveform and is explored in the subsequent section. 
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A. Fourth post-Newtonian dynamics 

In Sec. IIV XI we introduced the parameter 65 which encapsulates the unknown physical effects in the dynamics of 
the two bodies at orders higher than the 3PN order. The most relevant quantity that it affects is the potential A(u; 65 ) 
which occurs in the effective one-body Hamiltonian and the effective metric. Among other things A(u) governs the 
rate of inspiral of the bodies and therefore the phase of the waveform. 

We begin our visual comparison by plotting the effective potential A(u) at various PN orders including the 4PN 
order for two extreme values of 65 (b 5 = —50 and 65 = 500). Recall that u = 1/r ~ GM/ |xi — X 2 I and therefore 
u —> 0 denotes the region when the two bodies are infinitely separated and u ~ 0.5 denotes the region when the 
two black holes are “touching” each other. The sensitivity of ground-based interferometers is best in the frequency 
range 40-400 Hz. For a candidate system of total mass M = 20M 0 this frequency range corresponds to a range for 
the frequency-related variable v, Eq. 14.2611 , that begins at v = 0.2313 and terminates beyond 'cf s I ( ’ N = 0.4457 [ill or 
/iso = 286 Hz. Equivalently, this implies a range for the radius-related variable u [using Eqs. 13.1211 13.1311 to connect 
v and u], starting at u = 0.05335 or r = 1/u = 18.7. The gray-shaded region corresponds, when the total mass 
M = 20M 0 , to the frequency band [65,235] Hz, centered at 150 Hz, in which the SNR accumulated for inspirals is 
more than 80% of the total SNR in the entire LIGO band. For the system (1OM 0 ,1OM 0 ) the above frequency band 
corresponds to 0.2719 < v < 0.4174, equivalently 0.07372 < u < 0.1777 (or 13.56 > r > 5.63). The dashed vertical 
line at mi so = 0.2065 near the shaded region corresponds to the radial coordinate ri so = 4.84 at which the system 
reaches the last stable circular orbit. 

It is important to note that the (dashed) vertical line at ui so is invariant for systems of different masses but the 
shaded region will change with the total mass, moving to the right with increasing mass. The sensitivity of the 
instrument is best for those systems for which the LSO is close to (M ~ 20M 0 ) or within (M ~ 30M 0 ) the shaded 
region. 

From Fig. ^ we draw three important conclusions: (a) The potentials predicted by the two extreme values of 65 
used in our study encompass the variations implied by the second and third post-Newtonian orders, (b) In the region 
where the detector is most sensitive to binary black holes the agreement between the different models is pretty good, 
(c) Even consideration of extremely large positive values of the 4PN parameter 65 has little effect on the function 
A{u). [This is due to the fact that, after Padeing, the function A(u) has a limit when 65 —> 00 .] Even variations 
beyond reasonable values at 65 = —50 and 65 = 500 lead to an effective potential that is within the range of variation 
caused by different post-Newtonian orders. These observations already indicate that we should expect the fiducial 
EOB template to mimic flexed waveforms reasonably well. 


B. Unknown third post-Newtonian energy flux 9 


As mentioned in Sec. HTIthe gravitational energy flux at third post-Newtonian order has one undetermined parameter 
9. We also argued in Sec. ITVTI that the magnitude of 6 should be of order 1. Since the flux plays a crucial role in 
the phasing of the waves it is important to measure the effect of this parameter. In other words what fraction 
of the signal-to-noise ratio will be lost by setting 6 = 0 in our templates while in reality 9 is different from zero? 
Obviously, the answer depends on the extent by which 9 is different from zero. We vary 9 from —5 to +10 and plot 
the Newton-normalized flux as a function of the invariant velocity parameter v. As in Fig. ^ here too the shaded 
region corresponds to a frequency band [65, 235] Hz around 150 Hz corresponding to the range 0.2719 < v < 0.4174 
when the total mass M = 20M 0 (as assumed in Fig.^J. The dashed vertical line is of course the “velocity” at the last 
stable orbit i>i so = 0.446 corresponding to two systems of equal mass (independently of the value of the total mass). 

Figure [2 indicates the extent of variation caused by changing the value of 9. Clearly, negative values of 9 have a 
larger impact on the flux than positive values. Indeed, 9 = —5 leads to much greater variation in the flux than even 
9 = +10. The main message from Fig. Q is that by varying 9 over the range —5 < 9 < 10 in our study of faithfulness 
and effectualness we would in effect take into account the possibility that the real gravitational wave flux be rather 
different from that assumed in the effective one-body approximation. (Note, however, that the differences within the 
most relevant shaded region are only ~ ±10%.) The variation in 9 we consider is far greater than the variation of ±2 
in the parameter 9 considered in Ref. 0. Note that 9 is related to our 9 via 


Indeed our range corresponds to —3.49 < 9 < 11.5. 


1320 


(5.1) 






17 



u=1/r 


FIG. 1: The potential A(u) is plotted as a function u = 1/r ~ GM/|xi — X 2 I at various PN orders. By varying the 4PN 
parameter bs we more than cover the behaviour of both the 2nd and 3rd post-Newtonian orders. In all the Figures, for a total 
mass M = 20 Mg, the gray-shaded region corresponds to the frequency band [65,235] Hz, centered at 150 Hz, in which the 
signal-to-noise ratio accumulated for inspirals is more than 80% of the total SNR in the entire LIGO band. The corresponding 
range in u is 0.07372 < u < 0.1777 and r is 13.56 > r > 5.63. The dashed vertical line at ui so = 0.2065 near the shaded region 
corresponds to the radial coordinate ri ao = 4.84 at which the system reaches the last stable circular orbit. 


C. Flexibility parameter cp 

In Fig. 0 we plot the standard P-approximant flux at 3PN order in the equal mass case (i.e., 77 = 0.25, solid line). 
In the notation introduced in Sec. CVBl this curve corresponds to cp = 0. The effect of changing the location of 
the pole is shown by plotting the energy flux for four values of the parameter cp. We note that the curves change 
monotonically as the value of cp is changed, moving to the right for negative values of cp and to the left for positive 
values. We have changed the location of the pole by rougly 50% on either side of its nominal value. As mentioned 
before this amounts to varying the light ring value of v in the range 0.4605 < ui; g ht ring < 1.3814. No calculation we 
are aware of suggests a larger variation in the location of the pole than considered here. The curves do show a rather 
large variety indicating that the location of the pole is a priori as important as the other two parameters discussed 
before. 


VI. THE SPAN OF THE 3PN EOB BANK: OVERLAP OF FIDUCIAL TEMPLATE WITH FLEXED 

WAVEFORM 


A. Faithfulness and effectualness 


The ultimate tool for testing the robustness is of course the overlap of template waveforms with flexed waveforms. 
Given a fiducial template T(t]pk ) and a flexed signal S(t;pi,n a ) their overlap O is defined as 


0(T, S) = 


(T,S) 


y/{T,T)(S, S) 

where the scalar product is defined as usual by the Wiener formula 

df 


(X,Y) = 2 


'o S h (f) L 


X(f)Y*(f)+X*(f)Y(f) 


( 6 . 1 ) 


( 6 . 2 ) 
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r|=0.25 



FIG. 2: Variation in the (Newton-normalized) energy flux emitted by the system due to the 3PN parameter 6 being different 
from zero. Clearly, negative values of 6 have a greater effect on the behaviour of the flux as compared to the positive values. 


Here a denotes the Fourier transform of function a(t), that is, a(f) = f^° aft) exp(— 2mft)dt, a*(f) denotes complex 
conjugation of a(f) and Sh(f) is the (one-sided) noise spectral density of the detector. In computing overlaps we use 
the initial LIGO noise spectral density of Ref. [ 3 ] given by: 


S h {f) = 1.44 x IQ” 46 


( 4 . 64 /) ~ 56 



- 4.52 

+ 3.25 + 2 



Hz’ 


(6.3) 


where /*. = 150 Hz. Since the noise curve rises very steeply at low frequencies the lower limit of the integral in Eq. 
(IQ) does not have to be zero. It suffices to choose a lower limit of 40 Hz so as not to lose more than 1% of the 
overlap for binaries with total mass M ~ 50 Mq. 

Faithfulness T is defined as the overlap maximized only over the extrinsic parameters of the template, which in our 
case are simply a reference time t re f at which the template waveform reaches a certain frequency (say 40 Hz) and the 
phase (f >ref of the signal at that time: 


T = max \0(T,S)\. 

^•ref ?0ref 


(6.4) 


Effectualness £ is defined as the overlap maximized over not only the extrinsic parameters but all the intrinsic 
parameters as well, which in our case are the two masses mi and m 2 of the binary: 


£ = max \0(T, S)|. 

^ref j^ref j^T-2 


(6.5) 


B. Matched filtering and signal-to-noise ratio 

In searching for signals of known shape, such as chirping radiation from black hole binaries, one employs the method 
of matched filtering. For signals of known pattern, matched filtering is, in Gaussian noise background, a statistically 
optimum strategy in which a data analyst computes the cross-correlation of the template waveform T(t) with the 
detector output X(t). The analyst will not know before hand when the signal arrives or what its parameters are. 
Therefore, it is necessary to take several copies of the template corresponding to different parameter values pk and 
compute the correlation of each of those templates with the detector output at different time-steps to,ti,... ,tk, ■ ■ ■ ■ 
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r|=0.25 



FIG. 3: The (Newton-normalized) energy flux is plotted for different locations of the pole parametrised by cp. We vary the 
location of the pole by about 50% on either side of its nominal value predicted by the second post-Newtonian binding energy. 
The range of cp is perhaps far greater than what one could expect on physical grounds and causes a great variation in the flux 
function. Note that for cp = 0.5 the pole is moved to v po ie = 0.4605 which is near but still beyond the LSO. 


If detector output contains a sufficiently strong signal resembling one of the template waveforms then the cross¬ 
correlation will exceed the rms value of the correlation by a large amount, thereby generating a trigger for the analyst. 
It is well-known that the signal-to-noise ratio p of a template T with the detector output that contains a signal S of 
known shape is given by 


p=^ = -^0t = \O(T,S)\(S,S) 1/2 . (6.6) 

Here (S, S') 1 / 2 is the signal-to-noise ratio which would be achievable by matched filtering if the detector output 
were correlated with the exact replica of the signal hidden in the noise modulo the amplitude which is irrelevant. 
The above equation tells us that when we do not know the exact shape of the signal, the signal-to-noise ratio gets 
degraded and only a fraction equal to the overlap of the template used in the search with the exact signal expected 
to be present in the detector output is what is recoverable. Thus, while matched filtering is an excellent technique 
in detecting signals buried in noise, dephasing of the template relative to the signal can quickly degrade the quality 
of the output. This can be readily seen from Eq. (IQ) where the Fourier amplitudes of the template and the signal 
are multiplied together before being integrated over the frequency. These amplitudes coherently add up only when 
the phase of the template coincides with that of the signal at all points in the frequency space. Even a small initial 
difference in phase can accumulate and kill the integral since the signals last for a large number of cycles—more than 
60 cycles for (10Mq,10Mq) black hole binaries and up to 1000 for NS-NS binaries in the sensitivity band of the 
LIGO interferometers. This is the motivation for building template waveforms that are as close to the true general 
relativistic signal as possible. Note, however, that the weighting of the overlap integral Eq. itiPll by the inverse of 
the noise spectral density means that not all cycles in the signal and the template, are equally important. Ref. [3l| 
introduced the concept of “useful cycles” as a measure of the effective number of cycles which dominate the overlap 
integral. For instance in the case of a (10Mq,10Mq) system the number of useful cycles is ~ 7.6. This means 
that it is crucial to model the phasing of the signal during the ~ 8 cycles around the peak of the SNR logarithmic 
frequency-distribution i.e. around fp ~ 165Hz, but that a less accurate model of phasing for the other cycles may be 
acceptable. 
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TABLE II: Robustness of 3PN EOB model with respect to parameters 65 , 9, £2 and cp. We give the faithfulness T, effectualness 
£, and the total mass M' and symmetric-mass ratio rf that maximize the overlap while finding effectualness. We vary each 
parameter in the nominal range but also quote extreme values of the parameters up to which an overlap of 95% or greater 
is obtained. Some parameters, for example 65 will be bounded either from below or above since the physical quantity it 
participates in may be irregular for such values thereby affecting gravitational-wave phasing. 


System —> 


(10M o , IOMq) 


(15M O ,15M 0 ) 

parameter .[ 

T 

£ 


T 

£ 

(M',V) 

6 s = -50 

0.9754 

0.9944 

(20.25, 0.2479) 

0.9728 

0.9953 

(30.41, 0.2494) 

65 = +50 

0.9431 

0.9850 

(19.76, 0.2499) 

0.9448 

0.9943 

(29.25, 0.2499) 

65 = +190 

0.8566 

0.9663 

(19.61, 0.2500) 

0.8855 

0.9829 

(28.92, 0.2500) 

65 = +250 

0.8359 

0.9633 

(19.63, 0.2485) 

0.8476 

0.9766 

(28.54, 0.2500) 

65 = +300 

0.8014 

0.9558 

(19.48, 0.2500) 

0.8363 

0.9753 

(28.46, 0.2500) 

65 = +500 

0.7849 

0.9498 

(19.46, 0.2499) 

0.8181 

0.9694 

(28.46, 0.2499) 

cp = — 0.2 

0.9988 

0.9999 

(19.97, 0.2500) 

0.9993 

0.9998 

(29.96, 0.2500) 

cp = + 0.2 

0.9941 

0.9989 

(20.17, 0.2469) 

0.9975 

0.9999 

(30.31, 0.2469) 

e = -5 

0.9358 

0.9949 

(20.27, 0.2500) 

0.9642 

0.9971 

(30.77, 0.2497) 

e = +10 

0.9946 

0.9982 

(19.97, 0.2500) 

0.9964 

0.9999 

(29.86, 0.2498) 

C2 = -2 

0.9999 

1.0000 

(20.01, 0.2497) 

0.9998 

0.9998 

(30.02, 0.2493) 

C 2 = +2 

0.9999 

1.0000 

(20.01, 0.2500) 

0.9998 

0.9999 

(30.08, 0.2490) 


C. Span of a template bank 

Stringent as it may sound, matched filtering is not totally restrictive when dealing with a template bank rather 
than a single template. To explain how a template bank is not as restrictive as a single template we introduce the 
notion of span of a template bank. In the geometrical language of signal analysis Bsl ] templates can be thought of as 
vectors—one vector for every set of values of the parameters p k ■ If the signal depends on K parameters then the set 
of all vectors run over a Tv-dimensional manifold. Equation (l6~2l serves as a scalar product between different vectors 
and induces a natural metric gkm on the template-manifold with the parameters serving as natural coordinates: 

dT 

gkm = (Tk ,T m ), T k = -—. (6.7) 

apk 

Though each template is itself a vector in the vector space of all detector outputs, the set of all templates do not form 
a vector space. Therefore, when dealing with the problem of constructing a bank of templates one is really working 
with only a subspace of the vector space. Moreover, one does not work with the full template space either but, like 
in quantum mechanics, with the set of rays , i.e. the set of vectors modulo their lengths, which can be realized as the 
set of normalized vectors. In other words, we work on the unit sphere in the initial vector space. When considering 
a submanifold on this sphere which is not the intersection of a linear space with the unit sphere the metric g km gives 
only a local approximation to the vector product of the larger space, but it does not endow the finite-dimensional 
submanifold with the correct projection of the metric structure of the larger space. 

In our search for gravitational-wave signals we choose a grid of templates on the template-manifold. If the templates 
are an exact replica of the expected signal then the density of the grid points is so chosen that no signal vector on 
the manifold has an overlap with the “nearest” grid point smaller than a certain fraction called the minimal match 
MM [2fll l.‘l(|. typically chosen to be either MM = 0.965 or possibly MM = 0.95: 

max \0(T, S)| > MM. (6.8) 

Template bank 

The above inequality will be satisfied not only for signal vectors on the template manifold but also vectors that are 
off the manifold but close to it. In other words, the template bank obtains minimal match for all signals located in 
an infinite dimensional “slab” around the AT-dimensional template manifold but sufficiently close to it. This “slab” 
defines the span of the considered template bank, say S[T, MM]. Here, for simplicity we introduce one flexibility 
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parameter at a time and explore successively the “slab” along the directions defined by each extra parameter. It 
therefore suffices to consider only those signals that live in a (K + l)-dimensional space around the template manifold, 
which particular K + 1, depending on the flexibility parameter in question. The span S[T, MM] of a template bank 
along a given flexibility direction is then defined as the maximum domain in the corresponding (K + l)-dimensional 
space within which the template bank T obtains a given minimal match MM. In this work we estimate this domain 
by computing the range of the flexibility parameters within which the minimal match is achieved between the fiducial 
template and flexed waveforms. 

When the template is not a true representation of the signal, the signal vectors run over a manifold that does not 
exactly coincide with the template manifold. What is required for signal detection is that the span of the template 
bank includes the signal-manifold. Of course, if the minimal match is sufficiently small then any template bank would 
span the signal-manifold. Successful signal detection, without undue loss of signals, requires the signal-manifold to 
be a submanifold of S[T, 0.95] (i.e., 95% minimal match) or, better, of S[T, 0.965] (i.e., 96.5% minimal match). 

Finally, let us note that while the span is defined with respect to a continuum of template bank in reality we will 
have to be content with a finite lattice of templates. Therefore, it is not guaranteed that the (maximum of the) overlap 
of a finite template bank with an arbitrary flexed signal within the span S[T , 0.965] will be greater than 0.965. If the 
template lattice is chosen such that the minimal match is at least 0.965 then the maximum overlap reached within 
the span S[T , 0.965] of the template bank might be reduced to 0.965 2 ~ 0.93. So the actual loss in the event rate 
might be 1 - (0.965 2 ) 3 ~ 20%. 


D. Span of third post-Newtonian EOB template bank 

The main conclusion of this study is that the standard 3PN EOB template bank without any additional parameter, 
spans the extensions separately implied by the seven flexibility parameters that account for the unmodelled effects in 
the EOB formalism affecting both the dynamics and radiation flux. This is demonstrated in Tables El and EH where 
we have considered two archetypal binaries expected to be observed by initial interferometers—the first consisting 
of a pair of 10M 0 black holes and the second consisting of two 15M 0 black holes. In Table UH we have explored 
the faithfulness and effectualness of the 3PN EOB template bank (our fiducial template) with respect to the four 
important flexibility parameters 65 , cp, 6 1 and (2 and in Table EH the same, but for the less important flexibility 
parameters /NonAdiab, /NonCirc, and /transition- 

First, we discuss the results obtained by varying the parameters over the range in which they are expected to 
lie. Clearly, the 4PN parameter 65 has the strongest influence followed by the parameters cp,9 and £ 2 . Indeed, the 
faithfulness is not always larger than the fiducial minimal match of 0.965 when these parameters take values over the 
range in which they are expected to vary. However, for (15M 0 , 15M 0 ) systems the effectualness does easily meet the 
usual requirement S > 0.965 for all values of b 5 [including 65 > 500, in view of the fact, visible on Fig.E that A(u; 65 ) 
monotonically reaches a smooth limit as 65 —> + 00 ]. In the case of (1OM 0 ,1OM 0 ) systems the situation is a bit more 
involved: ( 1 ) when 65 < 200 , the usual requirements on effectualness is met, but ( 2 ) when 65 > 200 , the effectualness 
drops slightly below 0.965. However, for the “plausible” value 65 = 250 (see Sec. [E] above) the overlap is still larger 
than 0.963, and even for 65 = 500 (and probably for any larger 65 for the reason mentioned above) the overlap is still 
as large as 0.95. Note also that, in many cases, faithfulness is itself larger than the minimal match and effectualness 
is close to 1 . 

In view of this special sensitivity to 65 , we explored in detail the 65 -dependence of overlaps and found a simple 
modification of the standard 3PN EOB templates that allows for meeting the desired requirement 8 > 0.965 for 
all values of 65 . If one constructs a fiducial template bank by using as EOB potential the “ 65 -flexed” function 
% 50 (it) = % 4 pn(u; 6 5 = 50) [instead of A 3 pn(w)], we have found that it leads to effectualness larger than 0.965 in all 
cases [and in particular for the (1OM 0 ,1OM 0 ) system and —50 < 65 < 250]. Furthermore, as illustrated in Table ITVl 
the span of this new fiducial template bank now extends over all the values of 65 : — 50 < 65 < 2000. In the case of 
the less important parameters /NonAdiab, /NonCirc, and /transition, we observe that the faithfulness is itself larger than 
our minimal match except when /transition = 1 for the (15M 0 ,15M 0 ) system. 

As mentioned in Sec. E3 we have also explored robustness beyond the range in which the flexibility parameters 
are expected to lie and yet achieve the required effectualness. In Table E3 we summarize the range over which the 
different parameters are expected to vary together with the range over which they can be varied yet maintaining an 
effectualness of 0.965. This table shows that the span of the 3PN EOB bank of templates (or for that matter, the 
65 = 50-flexed 3PN one) in the other flexibility directions extend well beyond the expected plausible ranges. 
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TABLE III: Robustness of 3PN EOB model with respect to parameters /NonAdiab, /NonCirc, and /transition- We give the 
faithfulness T , effectualness £, and the total mass M' and symmetric-mass ratio rf that maximize the overlap while finding 
effectualness. We vary each parameter in the nominal range but also quote extreme values of the parameters up to which an 
overlap of 95% or greater is obtained. 


System —> 


(10 M 0 ,1OM 0 ) 


(15M 0 ,15M 0 ) 

Parameter j. 

T 

£ 

(M',r/) 

T 

£ 


/NonAdiab = 1-0 

0.9976 

1.0000 

(20.04, 0.2491) 

0.9967 

1.0000 

(30.11, 0.2486) 

/NonAdiab — “hl-0 

0.9979 

0.9999 

(20.00, 0.2500) 

0.9964 

0.9999 

(30.01, 0.2496) 

/NonCirc = 1*0 

0.9998 

0.9998 

(20.09, 0.2473) 

0.9997 

0.9999 

(30.02, 0.2500) 

/NonCirc = “hl.O 

0.9998 

0.9999 

(20.10, 0.2472) 

0.9996 

1.0000 

(30.04, 0.2500) 

/transition — 0.5 

1.0000 

1.0000 

(20.00, 0.2500) 

0.9967 

0.9994 

(30.03, 0.2498) 

/transition = 1.0 

0.9877 

0.9878 

(20.01, 0.2500) 

0.9531 

0.9606 

(30.26, 0.2498) 


E. How could the match be good when the flux functions look so very different? 

In Sec. El we noted that the behaviour of the energy flux Fp r (v) jF m (v) could be significantly different from their 
usual behaviour when the flexibility parameters are set to extreme values in their expected range. When the flux is 
so different how is it still possible to achieve good effectualness? 

The answer lies in several aspects of the problem: First, one should note that, after factorisation of the crucial 
“quadrupolar flux” oc v 10 , the changes in the Newton-normalized flux are less than 10%. Second, one should remember 
that for the massive binaries considered here, the number of “useful” gravitational-wave cycles |3l| corresponding 
roughly to shaded regions in Figs. [5] and [3] is quite moderate (< 8). Third, one should note that one of the crucial 
things affected by the flux is the total chirp time, or the duration, of the waveform. For instance, if the flux increases 
more rapidly when one of the flexibility parameters is nonzero, as in the case of 9 = —5, then the system loses energy 
more rapidly and therefore the waveform lasts shorter. However, this shortening of the waveform can also be achieved 
by making the binary heavier (or lengthened by making asymmetric binaries of the same total mass or simply lighter 
binaries). Recall that the effectualness is obtained by maximizing the overlap over both the extrinsic and intrinsic 
parameters. Thus, in the process of maximization one can absorb a change in time-scale by choosing a binary of 
different total mass and mass ratio. This explanation is borne out by a comparison of the trend of the curves in Figs. 
□ □ and El with the corresponding rows in Table GH For instance, effectualness for 9 > 0 (smaller flux at a given v 
than when 0 = 0) requires a system of total mass lighter than the original system while for cp > 0 (greater flux at a 
given v than when cp = 0) one requires a heavier system than the original one. 

A final comment: The insensitivity of the effectualness to the location of the pole can be interpreted to mean that 
the factorisation by the pole is not as crucial an element of the gravitational-wave flux resummation as perceived 
in E2- This suggests that it would be interesting to study the performance of templates which do not use such a 
factorisation of the flux function. We leave this study to future work. 


F. Systematic versus statistical errors in the estimation of parameters 

Together with the value of the effectualness, the tables also show the template parameters that obtain the maximized 
overlaps. A quick inspection reveals that the symmetric mass ratio rj of the template that obtained the maximum 
match is either equal to the actual value of 1/4 or when different from 1/4 the fractional difference is less than 
0.1%. This is probably explained by the following: We study templates as functions of mi and m 2 . But the function 
T(mi,m 2 ) is invariant under the permutation mi <-> m 2 and therefore the overlap 0(mi,m2) = (T(rn\, m 2 ), S) 
always reaches an extremum (along the lines mi + m 2 = constant) at mi = m 2 , i.e. at r] = 1/4. If these extrema are 
all maxima (as a function of the ratio mi/m 2 , for a fixed value of mi + 7712 ), the real maximum of the overlap must lie 
somewhere along the “ridge” mi = m 2 . Note, however, that there might as well be domains of parameter space where 
the overlap 0 (mi,m. 2 ) reaches a minimum along the ridge at m± = m 2 (at fixed mass scale mi + m 2 = constant). 
The total mass is different from the true total mass at worst by about 1.5%. These percentages, of course, do not 
give us a measure of the accuracy of estimation of the parameters, rather they tell us the extent of bias induced in 
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TABLE IV: The natural range for the flexibility parameters expected on physical grounds is shown together with the actual 
span for which the effectualness is greater than 0.965 for equal mass binaries of total mass 20 Mq and 30Mq. As discussed in 
the text the span in the 65 -direction refers to a bank of templates constructed with slightly modified EOB potential: Aso(u). 


Range/Span 

&5 

cp 

6 

C 2 

/NonAdiab 

/NonCirc 

/transition 

Expected Range 

[0, 250] 

[-0.2, 0.2] 

[-5,10] 

[-2,2] 

[-1,1] 

[- 1 , 1 ] 

[0,0.5] 

Span S(T, 0.965) 

[-50,2000] 

[-0.5, 0.5] 

[-5, >1000] 

[-5,100] 

[-20,15] 

[-50, 50] 

[ 0 , 1 ] 


the estimation. Since our template bank contains waveforms that are not exactly the same as the “true” signals, the 
parameters that maximize the overlap are different from the real values meaning there is a systematic error in the 
estimation of parameters and the percentages we have quoted are upper limits on the systematics. 

The statistical errors in the estimation of the intrinsic parameters m\ , m 2 are determined by the shape of the level 
contours of the overlap function between the template T(mi,m 2 ) and the signal, that we assume here to be part of 
the template bank: S = T(m f l ) , m®). For high signal-to-noise ratio this shape is an ellipse (the error ellipse) which is 
determined by the information matrix, i.e. the metric Eq. O- We have confirmed that even for the massive binary 
systems that we consider, these error contours are qualitatively well described by the analytical results of [H|, i.e. that, 
when represented in the (Ad, 77 )-plane, where Ad = is the chirp mass, they are highly elongated ellipses with 

major axis roughly along the 77 -direction and minor axis along the Ad-direction. When represented in the (7771,7772)- 
plane these error contours have the shape of thin crescents orthogonal to the diagonal mi = m 2 . A bad consequence 
of this fact is that there can be a large statistical error in the determination of 77 ]32), and therefore correspondingly 
large statistical errors in the determination of the individual masses. Only the combination Ad = rj 3 > 5 M might be 
reasonably free of statistical errors. A more detailed analysis is needed to assess the total errors combining systematic 
and statistical effects. 

We note in passing that a useful consequence of the highly elongated structure of the overlap contours is to allow 
a fast first-cut data analysis based on the lower-dimensional template bank defined by fixing 77 , e.g. to 1/4, and 
varying only M = mi + m 2 . The choice 77 = 1/4 would be sufficient to cover systems in a large domain of mass space 
around the diagonal mi = m 2 . It can then be complemented by considering a few other simple values of 77 . Each 
such lower-dimensional (approximate) template bank (corresponding to some value of 77 ) finally depends on only one 
universal function of one variable, the scaled phasing function (f>n(t), obtained by integrating, once for all, the EOB 
equations of motion expressed in terms of scaled variables t = t/M , r ~ |xi — X 2 I /GM, for a particular 77 , extending 
to EOB the idea of mother templates for the post-Newtonian model l 33 ll . The bank of templates is then built from 
the shifted and scaled function 4> v {{t — to)/M). This fact should simplify further the filter bank construction similar 
to the case of Newtonian signals which were expressible in terms of the universal function <f>N{i) oc — (—t) 5 / 8 . Indeed, 
the fact that this must be so is implicit in the nature of the template bank constructed by several authors (cf. second 
reference in [3(3 and Ref. [(3) who find that a small range of 77 is needed for a large range of the total mass. 


VII. CONCLUSIONS 

In this study we have explored the robustness of 3PN EOB templates. We introduced seven flexibility parameters 
that affect the two-body dynamics and radiation emission and varied each of them separately over a range that can 
be reasonably considered to be large enough to encompass unknown and unmodelled PN effects. The parameters 
introduced are: (a) a 4PN parameter 65 that alters the two-body effective metric and the EOB potential A(u). We 
conducted a special study of the structure of 4PN contributions to the Hamiltonian to estimate the plausible range 
of the parameter 65 measuring them (—50 < 65 < 250). (b) the unknown 3PN parameter 9 affecting the nature of the 
energy flux emitted by the system (—5 < 9 < 10). (c) a parameter (2 that changes 0(p 2 p 2 ) terms in the two-body 
Hamiltonian (—2 < C 2 < +2). (d) location of the pole in the energy flux controlled by a parameter cp. (e) and 
three parameters /NonAdiab, /NonCirc, and /transition, that are varied so as to inflict at least a factor two change in the 
modelling of non-adiabatic effects, non-circular effects and the transition from inspiral to plunge, respectively. 

We then compared the faithfulness and effectualness of standard fiducial EOB templates (that is, EOB templates 
in which all the above parameters are set equal to zero) with the flexed signals obtained by switching on the flexibility 
parameters one at a time. Based on the study conducted in this work we find that the third post-Newtonian EOB 
templates lead to effectualness larger than 96.3% in all cases (when 65 < 250). For the (10Mq, 10 Mq) systems, and 
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65 larger than 190 the effectualness drops below the usual requirement S > 0.965, though it remains very close to it, 
being larger than 0.96. (Even when b 5 gets very large the effectualness never drops below 0.95.) 

There are two ways of improving this situation linked to the special sensitivity to 65 . One way is to augument the 
standard 3PN bank of templates [based on A 3 pn(m)] by (when it is needed) a second bank of templates, based on 
Aioo(u) = A 4 pn(u ;&5 = 100). We have checked that this “doubled” bank of templates allows one to span all values 
of 65 with overlaps better than 0.985. A second way (which minimizes the total number of templates needed) is to 
work in all cases with only one specific 65 -flexed bank of templates [namely the one based on the “intermediate” EOB 
potential A$o(u) = A 4 pn(u ;65 = 50)]. The remarkable agreement between numerical and analytical descriptions of 
circular orbits near the LSO [ 3 , [nj, suggests that it might be possible soon to use numerical simulations to map in 
detail the EOB potential A(u) near the LSO. Hopefully this might lead to a numerical estimate of the value of 65 
thereby sharpening the preferred choice of bank of templates. 

There is a caveat in the current evaluation which we must bear in mind namely that we vary the flexibility parameters 
only one at a time. It is possible that the actual physical signal has more than one of the flexibility parameters non¬ 
zero. In that case our fiducial templates might not be able to achieve the desired span. This is because the ‘shifts’ 
in the fiducial template parameters needed to separately correct for the various non-zero flexibility parameters might 
be in different ‘directions’. An attempt to define a rather phenomenological signal was the main focus of Ref. 0. In 
this paper, we have focussed on a different (minimalistic) attitude and have not looked into these matters. We intend 
to address this issue in a future work. 

To conclude: the 3PN EOB templates (possibly suitably 65 -flexed) are good models to use for black hole binary 
searches in the interferometer gravitational-wave data because their “span” in signal space seems large enough to 
encompass most of the plausible modifications one can think of making in the current EOB framework. Moreover, as 
we emphasized, the highly elongated shape of overlap contours in the (M, ry)-plane suggests the interesting possibility 
to drastically reduce the number of EOB templates by using a small number of universal phasing functions (j> Vi ( t ) 
with a small discrete set of values of rj prominently including r) = 1/4. 
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